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Abstract. We address a possible relation between the expectation value of the Polyakov loop in pure 
gluodynamics and full QCD based on Polyakov Chiral Quark Models where constituent quarks and the 
Polyakov loop are coupled in a minimal way. To this end we use a center symmetry breaking Gaussian 
model for the Polyakov loop distribution which accurately reproduces gluodynamics data above the phase 
transition in terms of dimension 2 gluon condensate. The role played by the quantum and local nature of 
the Polyakov loop is emphasized. 
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The phase transition of QCD matter at finite temper- 
ature from hadrons to a quark-gluon plasma was estab- 

, lished long ago [DI21I3] (for a review see e.g. HJ). Tnc 
precise definition of the phase transition requires a proper 
identification of the relevant order parameters. In the heavy 
quark limit, the Polyakov loop vacuum expectation value 
signals the breaking of the center symmetry correspond- 
ing to the deconfincment phase transition when changing 
from zero to one jS] . In the limit of light quarks the chiral 
condensate determines the restoration of chiral symme- 
try when the chiral condensate vanishes above the critical 
temperature. In the real QCD case with dynamical mas- 
sive quarks both chiral and center symmetries are explic- 

1 itly broken, and neither the Polyakov loop nor the chiral 
condensate are truly order parameters, although a rather 
sharp crossover is expected across the phase transition for 
these quantities 0]. 

Polyakov-Chiral Quark Models allow to study the in- 
terplay between chiral symmetry restoration and center 
symmetry breaking [H] in a quantitative manner 7,8,9 
ITnillllll2llT^) . At zero temperature the constituent quark 
mass, M, is dynamically generated via the spontaneous 
breaking of chiral symmetry, inducing a exponentially small, 
~ e -M / T , breaking of the center symmetry at low tem- 
peratures |Ulll(Jllll| . This provides the rationale for keep- 
ing the Polyakov loop as an order parameter also in the 
unqucnched case. However, although the coupling of the 
Polyakov loop to quarks is rather unique, details regard- 
ing the postulated purely gluonic action differ 7,8,9,10 
II 111 12111. "S| . This additional information is beyond the chi- 
ral quark model capabilities and must always be postu- 
lated. In this regard, it seems natural to constrain the 
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gluonic action to reproduce pure gluodynamics results. In 
the present work we address a possible relation between 
the expectation value of the Polyakov loop in pure gluo- 
dynamics and full QCD based on the coupling to quarks 
generated by the fermion determinant within a Polyakov 
NJL model (PNJL). 

One of the problems one must also face when compar- 
ing models with lattice data has to do with the difficult but 
necessary renormalization of the Polyakov loop. Indeed, 
after renormalization the Polyakov loop falls outside the 
unitary group [TUTS] . The spontaneous breaking of the 
center symmetry above some critical temperature occurs 
already at the level of pure gluodynamics and is inter- 
preted as the signal of deconfinement |16j . Full dynamical 
QCD lattice simulations account, in addition, for an ex- 
plicit center symmetry breaking at low temperatures due 
to the presence of fermions j!7j . In a recent paper |18j we 
have shown that the behaviour of the Polyakov loop above 
the deconfinement phase transition can be naturally ac- 
commodated by a dimension-2 condensate, instead of the 
more conventional perturbative calculations. The quality 
of our fits leaves little doubt on the veracity of the descrip- 
tion of the lattice data. To our knowledge this is the first 
time that power corrections in temperature have been ad- 
vocated to describe QCD right above the deconfinement 
phase transition. This is in sharp contrast to the com- 
mon believe that logarithmic thermal corrections should 
be considered instead yielding to rather unnatural, in fact 
inaccurate, fits to the existing lattice results. 

Following our study on the role of dimension-2 conden- 
sate in the center symmetry breaking as a guideline |18| 
we make here a Gaussian ansatz for the purely gluonic 
contribution to the Polyakov loop distribution which re- 
produces remarkably well the lattice data for pure gluo- 
dynamics . Then, using the PNJL model [7051151111)11111 
I12U13| we may predict the value of the Polyakov loop in 
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the unquenched case and compare with full dynamical cal- 
culations. In our view this is the best possible scenario for 
a Polyakov-Chiral quark model to work. 

Let us explain how the dimension-2 gluon condensate 
can generate a power like behaviour of the Polyakov loop 
at temperatures slightly above the deconfinement phase 
transition in gluodynamics. In the Polyakov gauge the (ex- 
pectation value of the) Polyakov loop is defined as 



\N C 



L = ( ^-ti- r e igAo ^ T 



(1) 



where ( ) denotes vacuum expectation value and tr c is 
the (fundamental) color trace. Aq is the gluon field in 
the (Euclidean) time direction, Aq — J^^a^o.a, T a be- 
ing the Hermitian generators of the SU(iV c ) Lie algebra in 
the fundamental representation, with the standard nor- 
malization tv c (T a Tb) = 8 a i,/2. Note that in this gauge we 
have the invariance under the large gauge transformation 
gA -> gA a + 2ttT. 

If we expand the Polyakov loop at high temperatures 
we get the cumulant expansion 
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where odd terms have been discarded on the assumption 
that global colour symmetry breaking does not occur. Av- 
eraging over the N% — 1 colour gluonic degrees of freedom 
we get 



{Ao, ai ■ ■ ■ A),a 2fc ) 



((^o,a) 2fc )(iV c 2 -3)!! 
[N 2 + 2k-3)U 



(3) 



where S ai aah is the fully symmetrized product of the 
Kronecker delta symbols S a b- The SU(A r c ) Casimir invari- 
ant is C 2 = Ea T aT a = {N* - 1)/(2JV c )1jv . In the large 
N c limit only the contractions of adjacent indices (includ- 
ing cyclic permutations) survive. For instance, if we twist 
two adjacently contracted generators we get T a Ti,T a = 
—Tb/2N C so that (T a Tb) 2 is suppressed by two powers of 
N c as compared to T 2 T 2 . Then, neglecting these terms we 
obtain 
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Now, we expect ((^4o,a) 2fe ) to scale as N 2k , in which case 
and for large T the limit iV c — > oo with g 2 N c fixed is 
finite and nontrivial, at least in the absence of radiative 
corrections. Finally, if we assume vacuum saturation of 
condensates we can reduce further the expression since 



(A 2 k ) = (2k-iy.\(A^ + 



2\k 



(5) 



where n.v.c stands for non vacuum connected terms. After 
summing up the series we get the result 



exp 
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The exponentiation means that temperature is not neces- 
sarily very high since we summed up all the series and, 
actually, for small T yields L — ► 0. Of course, as one ap- 
proaches the true critical temperature one expects non 
analytical behaviour not encoded in the high temperature 
expansion 1 . On top of Eq. © one must take radiative 
corrections into account yielding, in fact, L > 1 (141115) 
because of renormalization effects at high temperatures. 
This is equivalent to multiply Eq. © by a slowly varying 
temperature dependent factor |18j . 

The previous behaviour suggests a Gaussian ansatz for 
the Polyakov loop distribution, i.e. the probability of hav- 
ing a configuration with a given value of (2. Unfortunately, 
not much is known about this object in QCD (see how- 
ever 51]). In the Polyakov gauge, where Ao is both static 
and diagonal we can write the piece of the vacuum wave 
functional depending on the temporal gluon fields as, 



Po (A 3l A s )=Ce-^ +A ^/ a2 



(7) 



where we assume that the gluon variables take any real 
values on the real axis. The normalization constant C and 
the width of the distribution a can be determined from 
the definition of the dimension 2 condensate 



(Al + A 2 8 ) = a 2 



(8) 



Although the previous Gaussian reproduces the result of 
the cumulant expansion of the Polyakov loop, it does not 
incorporate the known gauge invariance which corresponds 
to a periodicity condition on the gluon field distribution 
and the integration measure. 

Once the Polyakov loop distribution has been deter- 
mined for pure gluodynamics we set out to describe its 
expectation value when dynamical quarks are included. 
To do so we use the PNJL model. For a constant value of 
the Polyakov loop, S7, and the scalar field S — M (which 
we identify with the constituent quark mass) the quark 
effective action is given by 



T 
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(9) 



where we have only retained the vacuum contribution, so 
there is no contribution of meson fields. V is three dimen- 
sional volume and = +y / k 2 + M 2 is the energy of a 



1 Nonetheless, although Eq. © does not predict a phase 
transition one could still make an estimate of the critical 
temperature related to the point where L presents no curva- 
ture, i.e. the inflexion point which yields, T c 2 = g 2 (Aq) /6N C . 
On the other hand we expect the condensate to scale with 
N c and the number of Euclidean dimensions as g 2 (A 2 l }^ ~ 
4na(fj,)(N 2 — 1)£Mq Cd . This shows that the large N c limit 
of L is well defined for fixed T. Taking /Iqcd = 240 MeV 
and the scale /i = 2 GeV with the conversion factor fi = 
2ttT and iV c = 3 we get g 2 (Al) = (2.2 GeV) 2 in rough 
agreement with the determinations from lattice data. Then 
T 2 ~ 47ra(27rT c )(iV c 2 - l)^Q C D/6iV c yielding the estimate 
T c = I.LIqcd- 
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constituent quark with mass M. We define the Polyakov- 
loop averaged action 

e -r Q ( M ,T) = f d Q pQ ^ A& ) e -r Q {M t n, T ) (1Q) 



The value of M is determined by minimization of Fq(M, T) 
with respect to M, OTq (M, T) jdM = 0, which is the gap 
equation and determines M at a given temperature T, de- 
noted as M* — M(T). This procedure allows to compute 
the integration in meson fields at the mean field level. In 
addition, the relation between the (single flavour) chiral 
quark condensate, (qq), and the constituent quark mass, 
reads 



2G S (qqY 



-(M* 



(11) 



Any observable is obtained by using the temperature de- 
pendent mass M* . For instance the expectation value of 
the Polyakov loop is computed as 



{ > Jdn Po (A 3 ,A 8 )e- r ^ M '^ T ) 



(12) 



The integral in df2 in the case N c = 3 and in the Polyakov 
gauge can be computed numerically. At the mean field 
level the Polyakov distribution becomes a delta function, 
but this picture is obviously only valid at high enough 
temperature, where quantum effects are negligible. At low 
temperatures the mean field approximation can lead to 
wrong results. For example, for the expectation value of 
the Polyakov loop in the adjoint representation one ob- 
tains — 1/(AT? — 1), instead of zero, as has been observed 
in lattice data J5| • The color group integration solves this 
problem . Finally, we have taken into account here the 
local effects in the Polyakov loop by consider a confine- 
ment correlation volume V a — 8nT 3 /a 3 , which follows 
from the simplest correlation of two Polyakov loops 



Y J d 3 x J dntx c n{x)tx c nKy) = j>j d 



3 xe -a\3!-y\/T _ 



( 13 ) 

Higher correlation functions are expected to have a smaller 
contribution 

The result for the chiral condensate and Polyakov loop 
expectation value is presented in Fig.f^ We consider y/a = 
425 MeV and m u = rrid = m q = 5.5 MeV. As we see the 
agreement is reasonable although the quark model pre- 
dicts a shift in the critical temperature of about 50 MeV. 
This discrepancy is compatible with the model uncertain- 
ties discussed in our previous work [llj where the explicit 
center symmetry breaking of the pure gluonic action was 
disregarded. It is also compatible with the large rescaling 
advocated in |12|. Basically, the shift in the critical tem- 
perature reflects our lack of knowledge on the gluon action 
at intermediate temperatures, in particular, the feedback 
of quarks on the gluons. It would be most useful to develop 
a theoretical framework where this feedback is taken into 
account, as it is done, e.g. in the mean field approximation 
12 , while retaining the fundamentally quantum and local 
nature of the Polyakov loop rjTj • 
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Fig. 1. Temperature dependence of chiral condensate {qq}* 
and Polyakov loop expectation value L — (tr c O)/N c in the 
2-flavor PNJL model. We use the Gaussian distribution repro- 
ducing the Polyakov loop in the case of pure gluodynamics. 
Lattice data from |16| for gluodynamics (Nf = 0) and |17| for 
N f = 2 full QCD. Perturbative LO and NLO results for the 
Polyakov loop with Nf = 2 are shown for comparison. 
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